Dynamics of Domain Wall in a Biaxial Ferromagnet With Spin-torque 



Xin Liu*, Xiong-Jun Liu and Mo-Lin Ge 

Theoretical Physics Division, Nankai Institute of Mathematics, 
Nankai University, Tianjin 300071, P.R.China 

The dynamics of the domain wall (DW) in a biaxial ferromagnet interacting with a spin-polarized 
current are described by sine-gordon (SG) equation coupled with Gilbert damping term in this 
paper. Within our frame-work of this model, we obtain a threshold of the current in the motion of a 
single DW with the perturbation theory on kink soliton solution to the corresponding ferromagnetic 
system, and the threshold is shown to be dependent on the Gilbert damping term. Also, the motion 
properties of the DW are discussed for the zero- and nonzero-damping cases, which shows that our 
theory to describe the dynamics of the DW are self-consistent. 
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During the past decade or so, rapid advances have been witnessed in both theoretical and experi- 
mental aspects toward probing the novel mechanism of spin transfer torque, especially focused on the 
motion of the domain wall (DW) induced by the spin-polarized current. 

Current-driven motion of a DW was studied in a series works by Slonczewski Q| and Berger QQQ. 
In particular, the phenomenon that the electric current exerts a force on the DW via the exchange 
coupling was argued in 1984 2], and a spin-polarized current exerts a torque on the wall magnetization 
and also the motion due to the pulsed spin-polarized current are studied in 1992 3j, etc.. In the 
following years these interesting phenomena was observed by many experiments^ |(], 3. At the same 
time there are many theoretical efforts 0, 0, 0, 0, Q| to understand the microscopic origins. 
The theoretical works, however, still seem not satisfactorily to explain these phenomena and the 
intrinsical reason is still unclear.,. Ob the other hand, al, experiment* Qfl0 indicate a thread 
for the spin-current which induces the motion of the DW, i.e. the motion doesn't happen until the 
spin-current strength is lager than the value of this threshold. So, a central theory of the induced 
domain-wall-motion is to make clear the essential of the threshold, which we proceed to study in this 
paper. 

Very recently, the single-domain-wall as well as its dynamical properties has been observed directly 
|l4| , and another important phenomenon that current inducing a single DW switching exists in ferro- 
magnetic semiconductor has also been demonstrated in experiment |l5| . All these experiments offer 
more information for us to further understand it. At the same time, many discussions have been 
presented in the refs Q, 0, 0, Q| , especially the kind of spin-torque 

was proposed firstly in 16] and discussed in detail in the refs. 

HQ 

, where M is the ferromagnetic 

magnetization, M s is the saturation magnetization, P is the spin polarization of the current, [Ib is the 
Bohr magneton and j e is the electric current density. In this paper, considering the above type of spin 
torque, we study the dynamical properties of the DW in biaxial ferromagnet interacting with a spin- 
polarized current, focusing on the threshold of the current in the motion of the DW. The threshold is 
found to be dependent on the Gilbert damping term from the stability condition of the DW's static 
state, and the motion of the DW is shown to be a conclusion with the minimal energy of the system. 
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Now, let us consider a spin-polarized current propagating in a biaxial ferromagnetic material which 
has an easy x axis and a hard z axis. Here we take an example of a Neel wall and treat the spin 
configuration as uniform in y-z plane ^tj. We assume that the length of magnetization |M| is constant 
and the conducting electrons only interact with the local magnetization. In this case, the modified 
Landau-Lifschitz equation reads 

— = - 7 MxH e// + -Mx — +r (2) 

where 7 is the gyromagnetic ratio, 
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is effective magnetic field which includes the exchange field, anisotropy field and demagnetization field 
(note that there exists no external field when the current propagates in the ferromagnet, see, e.g. 
|l4,ll5j 'L H and H± are anisotropy fields which correspond to easy axis x and hard axis z respectively, 
and a is the Gilbert damping parameter. The spin-torque r given by Eq. Q only exists when the 
ferromagnet has a DW structure such as a microfabricated magnetic wire. The major object in present 
work is to study the threshold and the dynamics of the DW with Eq. (0) . To facilitate the further 
discussion, we derive the above equation with spherical coordinates (#,</>). By a straightforward 
calculation, we reach the following equations of and (f> 



sm6 - a \f<j>\ = /% 
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where 
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and 
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For convenience we note H± + AttM s as H± in the following derivation. To obtain the threshold, 
we firstly calculate the static solution to the Eq. Q and then examine its stability condition with 
perturbation theory. This is similar to the method used in the ref. |l6|| . But here, the magnetization 
M is a function of the space. 

For the biaxial ferromagnet, one has H± H, sin# ~ 1 and cos 8 ~ • Thus the static solution 
is governed by V$ = Ve = 0, i.e. 

bM+ 1C os6H ± = Q, (4) 
ox 

With the boundary condition </>(— 00) = and 0(oo) = 7T, one finds the form of the static kink soliton 
[ijj reads: 4> = 2 arctan[exp(±|)], where these solutions are referred to head-to-head DW (+ sign) 
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FIG. 1: (a) head-to-head DW, result of 0=2 arctan[exp(|)]. (b) tail-to-tail DW, result of cj>=2 arctan[exp( — |)]. 

and tail-to-tail DW (- sign) and £ = \J -jjjj — ^2^/ H± indicates that the spatial width of the soliton 
becomes narrower when the current strength increases (Fig.l). 

In the following part, we consider the head-to-head DW and proceed to calculate the threshold with 
perturbation theory. (This theory can also be used in a tail-to-tail DW and gets the same result.) For 
this we allow the static soliton solution develop a small dynamic factor, i.e. <p ~ > 4> '■ 

4>' = 2 arctan[exp((l + a(t))j + b(t))} (6) 

where b{t) is a small perturbation and a(t) is a higher order smallness of b(t) with a(t) ~ ( ) 2 (This 
can be derived from the dynamic solution of SG equation). Then, one has 5<j> = ft - <$> = ^b(t) and 
the equations of 8<f>, 89 can be written as 

d(5<j>\_ 1 (-bji£ jH x \(6<f>\_ 1 fi(6f\ (7] 
dt\S9J l + a 2 yabj^f -crfH± J \80 J 1 + a 2 \89 J { ' 

One can verify that the eigenvalues A of the matrix D is given by 

b j x 

Ai = 0, A2 = — r tanh(— ) — ajH± 

It is obvious that the stability condition of the solution <f> depends on the properties of the eigenvalues 
[2o| . Since Ai=0, the static solution can be steady only if A2 < at any position. For this one finds 
a critical value for the current 



_ ajH ± e I 2A 

P»b \Jm s H(1 + ^) { ) 

above which the static DW can not be stable any more, i.e. is the threshold of the current. Since 
the value of the function tanh x keeps approximate 1 from x — » 00 to the edge of the DW (we calculate 
this value with matlab, e.g. when x = 2£, one gets tanh(2) = 0.9640), the DW becomes unstable 
firstly at its edge. The pinning force |ld | is not considered in our derivation, thus not the pinning 
force, but the Gilbert damping term a is the exclusive factor for existing a threshold in our argument. 
Furthermore, if " S 1 - <C 1 (Because a may be very small, it's not contrary with the assumption that 
H± ;§> H), one finds the threshold is proportional to the H± and the wall-parameter \J as well 
Q 

Now, we proceed to reveal properties of the motion for the DW above the threshold. Also, we begin 
with Eq. (0. For the case H± H, sin 9 w 1 and cos9 w 0, the dynamic equation takes the following 
coupled form 

. B<j> 89 a<t> 

smt'— - a— = -bj- ^cos9H±_ (9) 

at ot ox 
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As the method used in |l8|, we differentiate equation © and substitute it into Eq. 1|10|) yields 

d 2 </> ( 2AfH± , 2 ,3 2 a 2 2 a^ 

^ - (— j£- - M^a + + 7 HH^T*** + a 7 H ± - = (11) 

From above formula one finds an additional term 2bj(j) xt (the third term of the left side) of Eq. 
lfTT|l is added due to the spin current. To facilitate the further discussions, we make transformations: 
t = tcos/3+nx sin/3, x — —tsin/3+KxcosP, tan 2/3 = bjn, where k = 1/ \J 2Al M f - — b 2 j = g^p^f ^C 
with ( > 1. Using these new coordinates, the Eq. l(TT)> recasts into 

' ; (S-S) + 7 2 ^^ + ^( CO ^-^ /3 g) = (12) 



cos 2/3 at' 2 ax' 2 ' dt' dx 

where ip = 2<j). To give a more intuitive discussion, we neglect the damping term in the following 

derivation. Then, under the boundary condition: oo) = ^>(oo) = 27r, we obtain kink soliton 
solution of the Eq. ijTTjl: 



2arctan[exp(\l^-^^-^^-(x - vt ))] 



2arctan[exp(t — c ° s ^ (nxtcos f3 — vsin/3) — t(sin/3 + v cos /?)))] (13) 

1 (1 — tr) 

where u is the velocity of the DW in x — t coordinate and a free parameter subject only to the 
restriction |u| < 1, and it is related to the velocity V by V= K ^ C o°g-"l1"n% • Witn tllis solution 
we can calculate the energy of the DW for the motion case. Deriving from the Lagrangian L = 
J d 3 x(M s <j>(cos9 — 1) — E), one can obtain the hamiltonian in the following form: 

H = Jd 3 xE = J d 3 x(2A((X79) 2 + sin(9(V</>) 2 ) - HM s (sm9 cosfi) 2 + H ± M S cos 2 9). (14) 

This result is valid no matter whether the DW is steady or not, and it indicates that there is no 
kinetic energy in the hamiltonian even for the propagating case of the domain wall. Substituting the 
Eq. ifKfll into lfTI)l. we obtain the DW energy E as 



E = 2 v / 2AM s Hcos 2/3({ 



cos (3 — v sin (3 vT~ 



Vl - v 2 cos2/3C 2 (cos/3- vsia/3) 
M s (bj(cosP — vsmfl) — (sin/3 + vcos(3)/k) 2 



-y 2 H ± 2A Vl-v 2 (cosf3-vsmP) 



} + C (15) 



where C is a constant independent on any parameters. Eq. (|15|l tells us that the total energy of 
the DW is the function of the parameters v and (3. The first term in the bracket of the right side 
of above equation corresponds to exchange energy, the second term corresponds to the anisotropy 
energy of the easy axis and the third one corresponds to the anisotropy energy of the hard axis. By 
a straightforward calculation we find that the total energy E reaches its minimum when v = tan (3. 
Therefore in this case, the DW structure can be described as 

cos 2/3 

<t>min = 2 arctan[exp( — - — (x - bjt))] (16) 
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FIG. 2: DW energy as a function of parameter v(\v\ < 1). The velocity of the DW V= K( "o° S /3-t,1° n % ■ Green 
square represents the energy of DW when V = and red hexagon represents the minimal energy of DW when 
V = bj(we take an example of (5 = -|). 



The above formula clearly shows that the distribution of 4>min is similar to the static soliton solution 
except for a velocity bj and M z — 0. So, if the damping term a = 0, the DW moves once the spin 
current is applied, and the threshold arrives at zero. These results are self-consistent with the former 
derivation that the threshold is proportional to the damping a (see eq. (0). 

Before sum up, we should point out that practically the velocity observed in experiment will be 
smaller than bj for the nonzero damping term. For a brief discussion on this effect one may assume 
the initial velocity of the DW is bj |l8| in the above-threshold case. With perturbation theory on the 
SG soliton, after Lorentz transformation to the soliton's rest frame z — t and noting that ip = 2</>, we 
obtain a series of equations similar to that in ref. [2l| : 

(f>(z,r) = (j) mm (z,T) +£(z,t), (17) 

here 4> m in (see eq. is described in z — r frame and \£(z, r)| <C 1 an d 

£(z,T) = ^ Ub (T)f b (z)+ ( dku k (T)f k (z), (18) 



7TT-irT+ sin ^ + iy (^-tan2/3-^ = 0, 19 

OT OZ i OT oz 

where fb(z) = 2sech(z), fk(z) = (2Tr)iuj^ 1 [k + itanh(z)] exp(ifcz), cu 2 =-f 2 HH± cos 2/3, uju = 1 + k 2 

3 

and v= a ~" H± ^° s 2 2,3 ■ Note that in Eq. i|19fl tan 2/3 is not a small value, we approximately obtain the 
amplitude of the translational mode u b satisfies 

du h {r) 



(It 



(8 tan 2/3(1 - exp(-i/r)) + exp(-^r) J /(r) exp(^r)dr) (20) 

with /(t) = i(7r/2)5 (tan2/3)i/ dk[kuk(T)]/ (w k si nh ^irk). This formula means that in z — r frame, 
the DW will get a velocity equal to (—dub/dr) |2lJ. From above result one easily find the velocity 
of the DW V should be smaller than bj. While we can not let t ~ oo to get the terminal velocity 
because for large r, this perturbation theory breaks down since Ub grows with r. In fact, theoretically 



G 



the motion will not stop as long as the spin current strength is larger than the threshold value, because 
the damping term will disappear once the velocity becomes zero and then according to present result 
this static solution is obviously a high energy state, which is not steady. 

The development herein is outlined as follows. Initially, no spin current is applied to the biaxial 
ferromagnet, the variable 6 — ir/2 and the DW can be described as static soliton solution of SG 
equation. Secondly, a current is applied while the DW still keeps static if the current strength is 
smaller than the threshold value, and its width becomes short in x-y plane. On the other hand, 
the magnetization M develops a component M z as is shown in Eq. (0}, the value of M z and the 
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FIG. 3: With spin-torque, the out-of-plane component M z as a function of ^-coordinate, is uniform in y-z 
plane. 

anisotropy energy of hard axis increase with time due to the increasing current strength till which 
arrives at the threshold (Fig. 3). Thirdly, for the above-threshold case, the static DW can't support 
the increasing of M z , i.e. this case will result in the propagation of the DW along the motion direction 
of the conducting electrons in the current. Then, neglecting the Gilbert damping term, we discuss the 
dynamical properties of Landau-Lifshitz equation by discovering a series of solutions which respectively 
correspond to different energies of the DW and find the velocity of the DW induced by the spin-current 
is bj, which corresponds to the minimal energy of the DW. Finally, for the nonzero-damping case, the 
velocity of the DW is found to be smaller than bj. With these derivations we show our result is self- 
consistent that the threshold depends on the Gilbert damping term of the coupled Landau-Lifschitz 
equation. 

In conclusion the dynamical properties of a single DW is studied in biaxial ferromagnet induced 
by a torque exerted via the spin-polarized current. With the perturbation theory on the static kink 
soliton solution to the ferromagnetic system, we obtain a threshold of the current in the motion 
of domain wall and the threshold is shown to be dependent on the Gilbert damping term of the 
coupled Landau-Lifschitz equation. Also, the motion properties of the DW are discussed in the zero- 
and nonzero-damping cases, which shows that our theory to describe the dynamics of the DW are 
self-consistent. 
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in Nankai University. It's these lectures that lead us to this attractive field. This work is supported 
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